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Abstract 


J  A  reliability  analysis  using  second- 
moment  approximations  is  conducted 
on  two  types  of  fire-exposed,  un¬ 
protected  wood  floors— conventional 
wood  joist  and  floor-truss  assemblies. 
A  methodology  is  illustrated  by  which 
the  probability  of  structural  failure  of 
a  wood  floor  assembly  can  be 
evaluated.  This  probability,  with  the 
probabilities  of  failure  of  other 
system  components,  can  be  used  in  a 
systematic  analysis  that  apparently  is 
a  viable  approach  to  realistic  analysis 
of  building  fire  safety. 

The  use  of  reliability  analysis  to 
compare  relative  fire  safety  of  dif¬ 
ferent  floor  components  is 
demonstrated.  A  procedure  for  in¬ 
troducing  new  components  into  the 
market,  based  on  a  concept  of  an 
equal  safety  index  calculated  for  a 
component  with  a  proven  inservice 
record,  is  discussed. 
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Introduction 

The  effectiveness  of  structural 
assemblies— walls,  floors,  ceiling 
roofs—  to  act  as  barriers  to  fire 
growth  is  measured  by  using  the 
American  Society  for  Testing  and 
Materials  (ASTM)  E-119  (3)3  fire  en¬ 
durance  test.  This  method  requires  a 
typical  loaded  assembly  be  exposed 
to  fire  and  the  time-to-failure  record¬ 
ed.  This  time  is  then  employed  to  rate 
the  assembly.  An  assembly  that  main¬ 
tains  its  integrity  for  30  minutes,  or 
more,  is  said  to  have  a  30-minute  fire- 
endurance  rating.  The  rating  pro¬ 
cedure  is  then  used  by  codes  to 
regulate  the  assembly  designs 
satisfactory  for  buildings  of  various 
occupancies.  The  Department  of 
Housing  and  Urban  Development 
(HUD),  for  example,  specifies  exterior 
walls  shall  have  a  20-minute  en¬ 
durance  rating  and  floors  a  10-minute 
rating  in  one-  and  two-family 
residences  (41).  This  is  believed  to 
protect  occupants  and  enable 
firefighters  to  safely  combat  dwelling 
fires. 

One  of  several  shortcomings  of  the 
present  fire-endurance  evaluation  and 
rating  system  is  that  it  does  not 
realistically  measure  the  fire  safety 


afforded  by  assemblies.  It  simply 
compares  the  fire  performance  of  one 
assembly  with  that  of  another  by  us¬ 
ing  a  single-performance  measure— a 
single  test  of  each  assembly  to 
measure  its  fire  endurance  under  a 
standard  fire  condition.  If  safety  is 
discussed,  soon  the  realization  is  of  a 
need  for  a  measure  of  the  risk 
associated  with  an  event  occurring. 
Hence,  the  present  fire-endurance 
rating  system  provides  little  informa¬ 
tion  on  the  safety  provided  by  struc¬ 
tural  assemblies.  This  lack  is 
somewhat  unfortunate  because  of  the 
potential  benefits  to  resident  owners, 
developers  of  new  assemblies,  and 
the  firefighting  and  code  enforcement 
community  that  could  be  realized  by 
having  this  information. 

Approach 

In  this  paper  a  risk-based 
methodology  is  presented  as  well  as 
its  application  for  assessing  the  fire- 
endurance  safety  of  two  unprotected 
light-frame  assemblies— the  conven¬ 
tional  joist  assembly  and  a  floor-truss 
assembly. 

In  a  floor,  there  are  three  primary 
modes  of  failure:  Structural  collapse, 
flame  penetration,  and  excessive 
temperature  rise  on  the  unexposed 


surface  (3).  In  a  system  analysis  the 
probability  of  occurrence  of  each  of 
these  events  would  be  required.  In 
one-  and  two-family  dwellings,  only 
criteria  of  structural  assembly  failure 
are  evaluated  (41). 

One  objective  of  this  research  was 
to  propose  and  illustrate  a  method¬ 
ology  by  which  the  probability  of 
structural  failure  of  wood  floor 
assemblies  could  be  evaluated. 

Another  objective,  perhaps  more 
important  for  the  immediate  future, 
was  a  method  by  which  new  com¬ 
ponents  could  be  introduced  into  the 
market.  The  underlying  premise  is 
that  a  new  component  can  be 
substituted  for  a  conventional  and 
code-acceptable  component  if  it  pro¬ 
vides  the  same  "degree  of  safety"  if 
exposed  to  the  natural  elements  such 
as  live  and  dead  load,  loads  due  to 
earthquakes,  fire,  and  so  on.  The  in¬ 
vestigations  reported  here  address 
the  fire-performance  aspects  of  floor 
assemblies. 
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Figure  1.— Family  of  curves  for  fire  endurance  of  wood  floors  (20) 

<r1'2  =  1/2(1  -  s)[1  -  (T I2\f  s)]J  where 

g  _  applied  load 

~  breaking  load  (extreme  fiber  stress  of  11,000  pounds  per  square  inch  at  maximum  load.) 
design  allowable  stress 
ultimate  stress 
s_=  DIB  (depth/breadth), 

T_=  t/20/A, 

t_=  fire  endurance  (min),  and 
=  cross-section  area  (in.2). 


Background 

One  rational  approach  to  ac¬ 
complishing  both  objectives  is  an 
analysis  using  probabilistic  engineer¬ 
ing  methods.  Probabilistic  engineer¬ 
ing  is  by  no  means  a  new  discipline. 
Ang  and  others,  as  of  1972  in  a  review 
of  literature,  reported  on  some  355 
research  papers  on  structural  reliabili¬ 
ty  (4).  Textbooks  have  been  published 
on  the  subject,  but  research  publica¬ 


tions  with  application  to  steel,  con¬ 
crete,  and  wood  engineering  number 
only  about  20.  The  use  of  this  same 
theory  in  fire  situations  has  been  sug 
gested  and  illustrated  by  researchers 
(5,9, 18,22). 

Fire  Endurance  of  Floor 
Assemblies— Deterministic 

In  reinforced  concrete  and  steel 
design  areas,  the  concept  of  fire 
design  engineering  as  opposed  to 


strictly  fire  tests  is  gaining  accep¬ 
tance  (1,2,12,34.40).  Wood  floor 
assemblies  are  qualified  or  fire  rated 
based  on  test  only,  and  these  data 
have  been  published  for  assemblies 
with  at  least  a  1-hour  tire  rating  (13). 
Various  researchers  have  been  pro¬ 
moting  a  combination  of  design  and 
testing  resembling  the  manner  in 
which  the  field  of  structural  engineer¬ 
ing  emerged. 

Sunley  (38)  reports  the  variability  of 
strength  of  timber  assemblies,  when 
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Figure  2.— An  idealized  exposed  floor 
joist  subjected  to  fire  on  three  sides. 
Subfloor  protects  top  side  of  joist. 
Although  it  is  known  bottom  corners 
round,  straight  boundaries  are  used 
as  an  approximation.  After  time, 
t^and  char  rate,  C,  depth  of  char 
equals  C  •  I  on  a  joist  of  height,  d, 
and  thickness,  h. 

(M  148  531) 


exposed  to  ASTM  E-119  conditions,  is 
less  than  the  variability  of  assembly 
strength  at  normal  temperature.  Addi¬ 
tional  support  for  this  is  evident  in 
the  findings  of  Schaffer  (33)  and 
Knudsen  and  Schniewind  (19)  on  the 
strength  of  small,  clear  specimens  of 
wood  at  elevated  temperature.  Sunley 
argues  further  that  fire  design 
engineering  should  be  no  more  dif¬ 
ficult  than  designing  for  other  load 
types. 

Equations  have  been  developed  to 
predict  the  fire  endurance  of  fire- 
exposed  joist  floor  (20,21).  The 
method  of  analysis  was  empirical, 
having  as  a  starting  point  the  flexure 
formula  tor  pure  bending.  In  each 
case,  the  solutions  to  these  equa¬ 
tions  are  obtained  by  graphical 
methods  (fig.  1).  These  equations  can¬ 
not  be  used  conveniently  in  a  pro¬ 
babilistic  analysis;  thus,  additional 
modeling  is  needed.  A  purely 
analytical  approach  to  predicting 
strength  at  elevated  temperatures 
would  be  most  desirable,  and  a  con¬ 
siderable  amount  of  basic  input  data 
exists  (34).  However,  this  type  of 


modeling  has  not  been  accomplished; 
thus,  simplified  strength  models  will 
be  used  in  the  analysis  that  follows. 

Analysis  of  Time-to- 
Failure 

The  first  step  of  a  probabilistic 
solution  to  an  engineering  problem  is 
to  identify  two  (preferably  indepen¬ 
dent)  random  variables,  one  of  which 
represents  a  load  effect  and  the  other 
a  resistance  effect.  These  variables 
are  normally  denoted  by  load  S_and 
resistance  R.  When  R^and  S^have  the 
same  units  and  the  probability  of  R^ 
being  less  than  S_can  be  interpreted 
physically  as  failure,  the  stage  is  set 
for  a  meaningful  solution.  In  a  fire 
situation,  if  fire  endurance  is 
associated  with  load  S^and  time-to- 
failure  of  the  component  is 
associated  with  resistance  R,  the 
stated  requirements  are  met.  Model¬ 
ing  the  time-to-failure  of  various  floor 
assemblies  and  the  fire  severity  to 
which  they  may  be  exposed  is  the 
first  required  step  to  further  estimate 
safety  of  assemblies. 


Fire  Duration 

The  prediction  of  fire  duration,  and 
more  generally  fire  severity,  is  itself  a 
research  topic  (6,11,23,29,30,39).  For 
this  analysis,  it  will  suffice  to  use  the 
approach  reviewed  by  Lie  (21).  For 
ventilation-controlled  fires,  in  which 
the  duration  would  be  the  longest,  the 
equation  that  relates  fire  duration  to 
available  ventilation  (e.g.,  window 
area  and  window  height)  is  given  by 
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where  W,  Ap,  A yy,  and  H  will  be  treated 
as  random  and  defined  as 
W  =  fuel  load  density  (kg/m*) 

Ap  =  floor  area  (mJ) 

Ayy  =  window  area  (m!) 

=  window  height  (m) 

The  constant  5.5  has  units 
kg  min  ''  m  '5/* 


Model  Building  for  an 
Exposed  Floor  Joist 

The  failure  during  fire  exposure  is 
assumed  to  be  caused  by  charring  of 
the  three  exposed  sides  of  a  joist; 
this  loss  of  section,  coupled  with  the 
strength-reducing  influence  of 
elevated  temperature,  causes  rupture 
of  the  joist.  Although  burn-through 
and  elevated  temperatures  of  the 
unexposed  surface  can  be  additional 
failure  criteria,  they  are  not  con¬ 
sidered  in  this  analysis.  (These  failure 
criteria  relate  directly  to  the  floor- 
subfloor  design  that  can  be  analyzed 
separately.)  Load  sharing  and  com¬ 
posite  action  are  not  accounted  for 
directly  in  the  analysis;  however,  they 
should  eventually  be  included  in  an 
experimental  verification  of  the 
model. 

A  typical  floor-joist  section  is 
shown  in  figure  2;  the  shaded  region 
shows  an  idealized  charred  area. 
Schaffer  (34)  reports  bottom  corners 
round  when  charring  occurs;  further¬ 
more,  the  radius  of  the  corners  can 
be  approximated  by  the  depth  of  char. 
To  account  for  this  rounding  by  the 
moment  of  inertia  would  complicate 
the  computations  in  the  analysis,  and 
it  is  clear  that  the  error  involved  by 
assuming  straight  boundaries  is  of 
minor  concern. 

By  use  of  the  flexure  formula,  an 
equation  can  be  written  to  quantify 
failure  in  a  fire  situation  as 


M  Y  (tf,C) 
Ktf.C) 


oB 


(2) 


where 

M  =  applied  moment  caused  by 
both  dead  and  live  loads 
(in.-lb) 

tf  =  time-to-failure  (min) 

Y(tf,C)  =  distance  to  extreme  fiber 

-  being  a  function  of  time-to- 

failure  and  char  rate  (in.) 

l(tf,C)  =  moment  of  inertia  about 

-  an  axis  midheight  the  remain¬ 
ing  uncharred  section  (in.*) 

o  =  an  exposed  joist  perfor¬ 

mance  factor  that  relates 
normal-temperature  strength 
to  high-temperature  strength 

B  =  joist  modulus  of  rupture  at 

room  temperature  (Ib/in.*) 
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This  model  is  similar  to  that  used 
by  others  for  large  beams  under  fire 
exposure  (20,21)\  it  also  neglects  any 
contribution  to  strength  by  the  floor¬ 
ing  itself. 

The  selection  of  the  model  needs 
some  justification.  At  room 
temperature  a  =  1,  and  the  model  is 
exact  with  the  thought  that  modulus 
of  rupture  is  an  idealized  linear  state. 
The  actual  state  of  stress  is 
nonlinear,  but  the  model  is  adequate 
for  design  especially  if  “depth  effect" 
is  considered.  In  a  fire  situation,  the 
nonlinearities  are  expected  to  worsen 
because  the  joist  cross  section  will 
not  maintain  a  uniform  temperature. 
The  movement  of  the  neutral  axis 
caused  by  a  nonuniform  modulus  of 
elasticity  (MOE)  resulting  from  a  non- 
uniform  temperature  profile  can  be 
expected.  At  the  same  time  e  com¬ 
pressive  and  the  tensile  strengms  of 
the  wood  fibers  are  reduced  because 
of  elevated  temperatures  [Schaffer 
(33)].  It  may  be  possible  to  model  the 
net  effect  of  the  mentioned  behavior 
if  coupled  with  the  normal- 
temperature  nonlinearities  of  bending, 
but  this  has  not  been  reported  in  the 
literature.  By  introducing  the  exposed 
joist  performance  factor,  a,  these 
unknowns  will  be  accounted  for  in  an 
empirical  sense.  In  addition,  the  fac¬ 
tor  a  will  account  for  the  rounding  of 
the  corners  and  to  some  degree,  load 
sharing  and  composite  action  of  the 
flooring-joist  assembly. 

By  referring  to  figure  2,  it  can  be 
seen  that  equation  (2)  can  be  re¬ 
written  as 

M  (d  -  C  tf)  12  =  oB 

(b  -  2C  tf)  (d  -C  tf)3  It 2  (3) 

where 

b_  =  initial  joist  width  (in.) 

dL=  initial  joist  depth  (in.) 


The  remaining  variables  have  been 
defined  with  equation  (2).  All  of  the  fac¬ 
tors  except  band  d_are  treated  as  ran¬ 
dom  variables.  By  rearranging  equa¬ 
tion  (3),  a  cubic  equation  In  time,  tf 
results  — 

—  =  bd2  -  2Cd  (d  +  b)  t, 

aB  +  C2  (b  +  4d)  tj  -  2C3  tsf  (4) 

Whereas  cubic  equations  can  be 
readily  solved  by  hand  calculation  or 
computer,  the  derivation  of  the 
statistics  of  the  variable  t^  would  be 
cumbersome.  This  led  the  authors  to 
investigate  the  error  introduced  in  tf 
by  dropping  the  cubic  term  and  simp¬ 
ly  solving  the  quadratic  equation  for 
tf.  Fortunately,  the  errors  introduced 
are  negligible,  ranging  from  1.58  per¬ 
cent  for  a  nominal  2  by  6  (38  by  140 
mm)  to  0.35  percent  for  a  nominal  2 
by  12  (38  by  286  mm).  It  must  be  em¬ 
phasized  that  this  approximation  was 
sho.vn  to  be  adequate  for  only 
nominal  2  by  6  (38  by  140  mm), 
nominal  2  by  8  (38  by  184  mm), 
nominal  2  by  10  (38  by  235  mm),  and 
nominal  2  by  12’s  (38  by 
286  mm). 

The  data  and  results  of  full-scale 
floor  section  (two  Douglas-fir  joists 
per  assembly)  tests  reported  by 
Lawson  (20)  were  used  to  study  the 
applicability  of  the  model.  A  question 
immediately  arises  regarding  the  ap¬ 
propriate  value  of  the  modulus  of  rup¬ 
ture,  B.  Because  the  test  for  the 
measurement  of  B^and  the  fire  test 
are  destructive,  it  is  impossible  to 
have  knowledge  of  both  properties  for 
a  single  piece.  At  this  impasse  the 
alternative  was  to  investigate  the 
mean  and  the  variability  of  the  pro¬ 
duct  of  the  two  variables  wand  B. 

This  was  done  in  the  following  man¬ 
ner: 

The  42  Douglas-flr  floor  assemblies 
tested  by  Lawson  were  divided  into 


Table  1.— Apparent  high-temperature  modulus  of  rupture,  aB,  for  each  grade  group 


Grade 

Sample 

size 

Apparent 

rupture 

strength 

Coefficient 
of  variation 
(COV) 

N 

*3’ 

°a  B 

Lbfln.’ 

kPa 

800 

20 

1,126 

7,764 

0.238 

1,200-1,600 

15 

1,126 

8,384 

.320 

1,800 

2 

1,264 

8,715 

442 

Clear 

5 

1,266 

8,729 

.199 

Combined 

42 

1,182 

8,150 

.271 

'  Constant  char  rate.  C.  ot  0  025  in./min  (0  0635  mm/mm)  was  used  in  calculation;  therefore.  S  values 
reported  are  an  upper  Bound  because  they  incfude  variability  ot  char  rate 


four  grade  groups— 800,  1,200-1,600, 
1,600,  and  clear— as  shown  in  table  1. 
The  data  shown  in  the  table  indicate 
little  difference  in  the  apparent  high- 
temperature  modulus  of  rupture  for 
the  four  different  grades.  More  sur¬ 
prising  are  the  calculated  coefficients 
of  variation  (COV).  From  data  col¬ 
lected  by  Hoyle  (15)  in  a  world  search 
of  lumber  data,  a  value  of  0.45  for  the 
COV  of  BJs  typical  for  Construction 
grade  joist.  With  all  grades  combined, 
the  Q^was  only  0.271,  which  shows 
E-119  exposure  (3)  apparently  does 
have  a  variance-reducing  influence  on 
BJust  as  Sunley  (38)  suggested. 

To  scrutinize  the  model  further,  an 
attempt  was  made  to  predict  the  fire 
endurance  of  two  Douglas-fir 
assemblies  tested  by  Son  (35).  Struc¬ 
tural  failures  occurred  at  11.63  and 
13.00  minutes  for  nominal  2  by  10  and 
2  by  8  joist  floors,  respectively.  By 
assuming  the  2  by  10  and  2  by  8  joist 
floors  of  similar  quality  as  those  of 
the  Lawson  report,  i.e.,  aB  =  1,165 
pounds  per  square  inch  (lb/in.2)  (8,032 
kPa),  the  predicted  time-to-failure  was 
2.58  minutes.  The  large  discrepancy 
between  actual  and  predicted  time  is 
attributable  to  the  difference  between 
the  live  load  levels  used  by  Lawson 
(20)  and  Son.  In  Son  s  tests,  the  joists 
were  stressed  to  100  percent  of  the 
allowable  design  stress,  as  is 
specified  by  ASTM  E-119,  whereas  the 
load  levels  used  in  the  Lawson  tests 
ranged  from  only  200  to  917  lb/in.2 
(1,378  to  6,322  kPa),  which  is  approx¬ 
imately  16  to  75  percent  of  the 
allowable  design  stress.  This  results 
in  failure  times  calculated  using 
figure  1  of  23  and  11.4  minutes, 
respectively,  for  a  nominal  2  by  8  joist 
floor  as  Son  tested.  Hence,  as  ex¬ 
pected,  the  lower  the  load  level,  the 
greater  is  the  time-to-failure;  thus, 
more  thermal  degradation  is  allowed 
to  occur.  Beams  more  heavily  loaded 
will  have  shorter  times-to-failure  and 
less  thermal  degradation.  The  right 
side  of  equation  (3),  representing  the 
strength  change  as  influenced  by 
temperature  rise,  cannot  account  for 
heat  accumulation  degrading  the 
cross-sectional  strength  unless  some 
measure  of  time  dependence  is  in¬ 
troduced.  By  making  several  data 
plots,  the  following  model  was 
developed  to  include  this  kind  of  time 
dependence  in  which  the  variables 
have  been  previously  defined. 


I 


4 


M(d  •  Ctf)  12 
(b  -  2Ctf)  (d  -  Ctf)3  /1 2 


B 


The  [(b  +  2d)  ytfj/(bd)  term  may  be 
viewed  as  a  time-dependent  geo¬ 
metric  factor  to  account  for  heat  flow¬ 
ing  into  the  cross  section,  bd,  through 
the  perimeter,  b  +  2d.  Although 
visual  inspection  of  the  data  plots 
suggests  using  this  model,  some 
variation  of  the  model  may  be  more 


suitable.  Again,  as  before,  the  cubic 
term  tf_of  equation  (4)  is  negligible.  A 
least-squares  nonlinear  regression 
analysis  was  conducted  on  five  varia¬ 
tions  of  equation  (4)  as  shown  in 
table  2.  The  models  were  fitted  to  the 
42  full-scale  floor  section  tests  of 
Lawson  (20). 

The  second  model,  a  single 
parameter  of  table  2,  is  seen  as  the 
preferred  model  because  the  residual 
standard  deviation  is  only  slightly 
greater  than  that  of  the  two- 
parameter  model  listed  first  in  the 
table.  The  range  of  predicted  times-to- 
failure  for  the  42  floor  assemblies 
was  7  to  29  minutes.  This  shows  the 


selected  model  is  a  reliable  predictor 
if  it  is  recognized  that  the  residual 
standard  deviation  of  2.57  minutes  in¬ 
cludes  the  variability  of  oB^as  shown 
in  table  1. 

Solving  equation  (5)  for  tf^by  omitting 
the  cubic,  results  in 

2Cd(d  +  b)  +  6MKy/B  (6) 
tf  =  - 

V{2Cd(d  +  b)  +  6MKy/B|  3 

•  4C!(b  +  4d)(bd3  -  6M/B) 

2C‘  (b  +  4d) 

where  =  (b  +  2d)/bd,  which  is  an 
explicit  expression  for  the  time-to- 
failure.  Time-to-failure,  tf,  is  compat¬ 
ible  to  the  previously  defined  fire 
duration  tg,  which  is  the  load 
variable. 

Results  predicted  by  equation  (5) 
were  compared  with  results  of  four 
floor  fire-endurance  tests  obtained  by 
the  NFPA  (National  Forest  Products 
Association)  (24,25,26,27).  The  actual 
times-to-failure  to  carry  load  versus 
those  predicted  are  given  in  table  3; 
the  predicted  times  are  consistently 
and  significantly  less  than  the  times- 
to-failure  of  the  whole  floor  assembly. 
This  deviation  can  be  explained  by 
the  model  parameter  derivation  based 
on  the  results  of  Lawson  (20).  Lawson 
conducted  fire-endurance  tests  of 
paired  Douglas-fir  joists  with  essen¬ 
tially  noncontinuous  floor  sheathing. 
The  NFPA  tests  are  of  assemblies  of 
many  joists  and  a  more  or  less  con¬ 
tinuous  floor  sheathing.  These 
assemblies  result  in  load-sharing  be¬ 
tween  joists  and  increased  load¬ 
carrying  capacity  of  the  sheathing 
that  is  not  similarly  reflected  in  the 


Table  2.— Least-squares  nonlinear  regression  analysis  on  five  different  variations  of  a 
time-to-structural-failure  model' 


Model3 

Parameter 

estimates 

Residual 

standard 

deviation 

Min 

M(d  -  Ctp/2 

B 

y0  =  0.296 

2.54 

(b  -  2CtfXd  -  Ctf)’/12 

y0  +  7,  Ktf 

y,  =  0.206 

_ 

B 

y  =  0.170 

2.57 

1  +  K  y  tf 

= 

B 

y0  =  120 

2.60 

1  +  KForitf 

y,  =  0.160 

B 

y  =  0.216 

2.74 

1  +  y  tf 

B 

y,  =  -0.324 

3.26 

1  +  K(y,  t(  +  y,  tf3) 

r,  =  0.0328 

'  In  each  variation,  42  data  points  were  used  to  estimate  either  one  or  two  parameters,  depending 
on  form  of  the  model  Parameters  must  carry  necessary  units  to  make  equations  dimensionally 
homogeneous. 

*  K  =  (b  +  2  •  dVtb  '  d). 


Table  3.— Predicted  and  actual  times-to-failure  for  NFPA  unprotected  floor  fire-endurance  tests  (24,25,26,27)' 


Sample 

Nominal 

size 

Applied 

joist 

moment 

Bll 

C 

Predicted 

•f 

Assembly— 

observed 

tf 

Joist3 

observed 

‘f 

In. 

In.-lb 

Lb/ln.3 

In./min 

Min 

Min 

Min 

No  2  Douglas-fir  "S-dry"  w/vinyl  tile, 
19/32  in.  plywood 

2x8 

19.054 

4,308 

0.0245 

5.42 

10.2 

5.0 

No.  2  Douglas-fir  "S-dry”  w/nylon 
carpet.  19/32  in  plywood 

2x8 

19.054 

4,308 

0.0245 

5.42 

12.86 

11.5 

No.  2  MG  southern  pine  "S-dry" 
w/vinyl  tile,  23/32  in.  plywood 

2  x  10 

31,017 

5,730 

0.03 

7.01 

13.34 

9.0 

No.  2  MG  southern  pine  "S-dry" 
w/nylon  carpet.  23/32  in.  plywood 

2  x  10 

31,017 

5,730 

0.03 

7.01 

12.06 

12.06 

'  NFPA.  National  Forest  Products  Association,  italicized  numbers  in  parentheses  refer  to  Literature  Cited  at  end  of  paper 
*  Douglas-fir  rupture  strength,  from  PS);  southern  pine  rupture  strength,  from  (7) 

1  Failure  time  tor  first  |Oists  to  fail,  not  assembly  failure  time 


5 


Table  4.— Fire  endurance  of  2  by  4’e  under  constant  tensile  load' 


paired  joist  tests  by  Lawson.  If  the 
time  is  noted  when  the  first  joist  rup¬ 
tures  in  the  NFPA  tests,  the  dif¬ 
ference  between  the  model-predicted 
results  and  those  actually  observed 
are  closer  in  two  of  the  tests  in  which 
such  failure  was  observed.  This, 
again,  illustrates  the  need  for  a 
degrade  parameter,  y,  or  other 
parameters  that  include  both  the  ef¬ 
fect  of  load-sharing  and  that  of  floor 
sheathing.  As  a  result,  these  types  of 
replicate  experiments  are  planned. 

Model  for  Exposed  Floor 
Truss 

The  lower  chord  of  a  floor  truss  is 
subjected  to  both  bending  and  ten¬ 
sion,  and  the  well-known  interaction 
equation  is  used  for  design  purposes. 


Douglas  fir,  coast 

Southern  pine 

Teat 

Failure 

time 

Test 

Failure 

time 

Min 

Min 

P  =  6,100  lb,  90  pet  Ft 

P  =  6,100  lb,  83  pet  F, 

1 

1  1 

11.20 

10.00 

2 

7.67 

2 

11.61 

3 

9.35 

3 

12.85 

4 

11.25 

4 

12  34 

5 

8.74 

5 

11.80 

Mean 

9.64 

Mean 

11.72 

Standard  deviation 

P  =  4,960  lb,  73  pet  F, 

6  1 

1. 57 

Standard  deviation 

1.08 

9.24 

7 

9.96 

8 

13.36 

9 

13.92 

Mean 

11.62 

Standard  deviation 

2.36 

fb  y.  .  '  Select  Structural  allowable  stresses  for  Douglasfir  and  southern  pine  from  1977  National  Design 

I  =  —  +  —  (7)  Specifications. 

Fb  Ft 


Here  fb  and  ^denote  applied 
stresses;  Fh  and  F(,  allowable  design 
stresses  in  bending  and  tension, 
respectively. 

As  in  previous  reliability  work  at 
the  Forest  Products  Laboratory  (37), 
this  interaction  equation  can  be 
modified  to  indicate  failure  (with 
some  reservations  discussed  in  the 
report).  However,  in  a  fire-exposure 
case,  one  parameter  needs  to  be 
estimated;  thus  some  slight  inac¬ 
curacy  in  the  neighborhood  of  the 
combined  stresses  associated  with  a 
floor  truss  will  be  corrected.  The 
failure  equation  for  fire  exposure 
would  read 


*=-+-  =  - : -  (8) 

B  T  1  +  g(b,d,tf,y) 

where  the  right  side  of  the  equation 
has  a  form  similar  to  that  for  the  ex¬ 
posed  floor  joist.  Function  ^accounts 
for  the  thermal  degrade  of  the  sec¬ 
tion,  and  its  arguments  are  later 
defined.  Bjs  the  modulus  of  rupture 
from  which  Fh  was  derived;  and  T^the 
ultimate  tensile  strength  property 
from  which  the  design  value  Ft  was 
derived. 

Because  four-sided  fire  exposure  of 
the  lower  chord  in  a  floor  truss  is 
critical,  expansion  of  the  interaction 
formula  for  this  case  is  given  in  the 
following  equation.  (It  is  assumed 
that  the  mode  of  failure  is  rupture  of 
the  lower  chord.  The  upper  chord  has 


only  three-side  exposure,  and  the 
webs  are  only  stressed  to  approx¬ 
imately  one-half  the  level  of  the  lower 
chord.) 

_ P _ 

(b  -  2C  tf)(d  -  2C  tf) 

T 


M(d  -  2V  if)  12 
(b  -  2C  tfXd  -  2C  tt)3  /12 
B 


1 

1  +  yK  t, 


(9) 


where 

K^=  2(b  +  d)/(bd) 

yj=  thermal  degrade  factor 

P^  =  axial  tensile  force  due  to  dead 
plus  live  load 

b  =  width 

d  =  depth 

=  char  rate 

t;  =  time-to-failure 

M  =  maximum  bending  moment 
caused  by  dead  plus  live  load 

B^  =  modulus  of  rupture 
T  =  ultimate  tensile  stress 


Analogous  to  the  floor-joist  case,  K  is 
the  ratio  of  the  lower  chord  perimeTer 
(or  surface  area  for  heat  transfer)  to 
the  cross-sectional  area  (or  volume 
for  heat  storage).  After  some 
manipulation,  a  cubic  equation 
results: 


!-8C3)t|  +  j 4CJ(b  +  2d)  +  —  ~  y-K  j-tj 
+  <j-2dC(d  +  2b)  -  -y-  yK  +  yC 


6MyK) 

u .  ^ 

-  ™Uo 

B  /  ’ 

1  T 

B  / 

(10) 


This  equation  will  later  be  solved  for 
tf  and  used  to  estimate  assembly 
reliability. 

Model  Parameters  for  an 
Exposed  Floor  Truss 

To  properly  estimate  model 
parameters,  test  data  on  2  by  4 
assembly  members  under  combined 
tension  and  bending  are  required.  Un¬ 
fortunately,  no  data  are  available  for 
this  purpose.  Schaffer  (31),  however, 
has  conducted  fire-exposure  tests  of 
constantly  tension-loaded  Select 
Structural  coast  Douglas  fir  and 
southern  pine  2  by  4  members.  The 
time-to-failure  was  recorded  (table  4). 


For  pure  tension,  the  failure  model 
reduces  to 


DOUGLAS -FIR  2«4 


P(i  +  r  K  tf) 

=  T(b  -  2C  tf)(d  -  2C  tf)  (1 

This  expression  is  easy  to  deal  with, 
because  it  is  only  a  quadratic  in 
After  solving  it  for  tf^the  results  are 

t  =  2T  C(b  +  d)  -t-  PyK 

V{  2T  C(b  +  d)  +  PyK 
-  16T  CJ  (Tbd  -P) 


Af /A|  *  (b  -  2Cmt)(d  -  2Cm*)/bd 


WITH  C_*  0.025  IN./MIN 


A,  /A,  =  ( b  -  2Cf  t)(d-  2Cf  t  )/bd 

WITH  Cf* 0.038  IN./MIN 


where  the  minus  root  is  the  mean¬ 
ingful  root. 

An  estimate  is  needed  for  the  mean 
tensile  strength,  T^  and  char  rate,  C, 
to  determine  yjor  the  available  ten¬ 
sion  fire  test  data^ 

An  estimate  of  TJs  available  for  in¬ 
land  Douglas-fir  Select  Structural 
[Hoyle  (14)].  Based  on  a  sample  size 
o'  30,  the  average  tensile  strength 
value  was  5,020  lb/in.J  (34,600  kPa) 
with  a  COV  of  0.388.  This  value  can 
be  compared  to  what  might  be 
calculated  using  normal  distribution 
theory  and  the  allowable  tensile 
stress.  From  the  1977  National 
Design  Specification  (28)  the 
allowable  tension  value  for  Select 
Structural  Douglas-fir/larch  is 
1,200  lb/in.'  Using  the  calculated  COV 
of  0.388,  the  mean  value  is  deter¬ 
mined  to  be  6,811  lb/in.J  by  the  follow¬ 
ing  formula: 


TIME,)  (MIN) 


SOUTHERN  PINE  2*4 


Af  /A|  =  (b-2Cmt)(d-2Cmt)/bd 

WITH  Cm*  0.030  IN. /MIN 


Aj/A^b-aCjtXd-ac^j/bd- 


WITH  Cf*  0.045  IN. /MIN 


A{/A.  =  100-6  981 


1,200  =  (T  -  1.645  *  COV  *  T)/2.1 


TIME,!  (MIN) 


This  value  is  significantly  larger  than 
the  actual  mean  of  5,020  that  illus¬ 
trates  the  non-normal  nature  of  ten¬ 
sion  data.  This  shows,  therefore,  the 
mean  value  obtained  from  lumber 
tests  should  be  used  whenever  possi¬ 
ble  for  a  grade  and  a  species  in  ques¬ 
tion. 

The  nine  Douglas-fir  fire  and  ten¬ 
sion  test  results  can  then  be  used  to 
estimate  ^for  Douglas-fir.  By  using  T^ 
=  5,020  and  C  =  0.0245  inch  per 
minute,  yjwas  estimated  to  be  0.113 
with  a  residual  standard  deviation  of 
the  time-to-failure  of  1.829  minutes. 


Figure  3  — Residual  cross-sectional  area  divided  by  initial  area  (pet)  of  coast 
Douglas-fir  (top)  and  southern  pine  (bottom)  nominal  2  by  4  (1-  5/8  by  3-  5/8) 
members  with  duration  of  exposure  to  ASTM  E-  119  fire  conditions.  Top 
curves  result  from  using  the  mean  species  charring  rate.  Cm,  of  0.025  and 
of  0.030  in. /min  for  large  sections  of  coast  Douglas-fir  and  southern  pine, 
respectively.  Bottom  curves  from  using  an  effective  charring  rate.  Cj  of  0.038 
and  of  0.45  in. /min  for  coast  Douglas-fir  and  southern  pine,  respectively. 
Straight  lines  are  linear  regression  fits  to  actual  data  in  which  lines  were 
forced  through  100  percent  at  time  equal  to  zero.  (Model  Y  =  px  was  fitted 
to  one  residual  areaNnitial  area.)  Equation  of  the  curves  is  given  by 
100(b  -2Ct)(d  -2Ct)/(bd). 
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Figure  4.— Floor-truss  design  subjected  to  test  conditions  of  ASTM  E-  119.  Up¬ 
per  chord  was  loaded  with  tanks  simulating  a  uniform  load  of  55 .1  Ib/lt1 
that  resulted  in  a  combined  live  and  dead  load  of  60  ib/ft2. 
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A  similar  estimate  of  the  fire- 
exposure  reduction  factor  £_can  be 
made  for  the  southern  pine  test 
results  of  table  4.  The  southern  pine 
lumber  tested  was  almost  clear  of 
defects  and  ungraded.  The  lumber  ap¬ 
peared  to  be  of  a  quality  at  least  as 
good  as  Select  Structural.  There  are 
data  [Hoyle  and  Maloney  (17)]  that 
show  high-quality  southern  pine 
(2400f  machine  stress  rated  (MSR))  is 
stronger  in  tension  than  is  high- 
quality  Douglas-fir  (2400f  MSR)— the 
ratio  in  strengths  being  1.24.  To  arrive 
at  a  mean  value  for  the  tensile 
strength  Tj>f  Select  Structural 
southern  pine,  the  5,020  lb/in.J  value 
for  Douglas-fir  was  multiplied  by  1.24 
to  yield  6,233  lb/in.J  for  southern  pine. 
Using  a  mean  char  rate  for  southern 
pine  of  0.03  inch  per  minute  (32)  in  a 
nonlinear  regression  analysis  of  the 
fire  test  results,  the  yield  is  a  yot 
0.0839  inch  per  minute  with  a  residual 
standard  deviation  of  the  time-to- 
failure  of  1.077  minutes. 

Both  species  values  of  the  reduc¬ 
tion  factor,  ^JO.113  for  Douglas-fir, 
0.0839  for  southern  pine)  are  substan¬ 
tially  lower  than  that  of  0.17  for  the 
floor-joist  assembly  described  in  the 
preceding  section.  The  reasons  for 
this  are  not  clear.  It  is  known, 
however,  that  smaller  sections  (2  by 
4,  compared  to  2  by  8  or  2  by  10’s)  are 
likely  to  undergo  more  rapid  reduction 
in  cross  section  than  are  larger.  The 
difference  for  2  by  4  sections  of 
Douglas-fir  and  southern  pine  are 
shown  in  figure  3  (top  and  bottom, 
respectively)  as  a  function  of  fire- 
exposure  time. 

The  developed  model  and  the  para¬ 
meters  may  be  used  to  estimate  the 
structural  failure  of  a  given  floor-truss 
assembly.  This  was  done  for  the  truss 
shown  in  figure  4;  the  lumber  of  the 
floor  truss  is  No.  1  Dense  KD 
southern  pine.  B  for  No.  1  Dense 
southern  pine  was  obtained  from 
table  2  of  Doyle  and  Markwardt  (7). 
Because  tensile  strength  data  for 
Dense  No.  1  were  not  available,  data 
for  No.  1  KD  southern  pine  were 
taken  from  Doyle  and  Markwardt  (8). 
The  published  value  of  5,706  lb/in.J 
was  adjusted  to  5,646  lb/in.2  to  reflect 
current  standards  by  Forest  Products 
Laboratory  personnel. 

The  truss  was  then  analyzed  with  a 
Purdue  Plane  Structures  Analyzer  (36), 
and  as  normally  done,  the  center 
panel  of  the  lower  chord  was  most 
highly  stressed  with  an  axial  force  of 
4,209  pounds  and  a  bending  moment 


Figure  5.— Relationship  of  applied  load 
assembly  of  figure  4. 
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of  140  inch-pounds.  A  summary  of  the 
input  parameters  can  now  be  given 
for  the  fire-tested  floor  truss. 

0.0839  inch  per  minute 
P^  =  4,209  pounds 
b^=  3.5  inches 
d_=  1.5  inches 

=  0.03  inch  per  minute 
M  =  140  inches-pounds 


to  time-to-failure  for  floor-truss 


B^  =  9,410  pounds  per  square  inch 
T_  =  5,646  pounds  per  square  inch 

Substitution  into  and  solution  of  the 
failure-model  equation  (10)  results  in 
a  time-to-failure  estimate  of 
11.2  minutes;  the  actual  failure  time 
resulting  in  a  fire  endurance  test  was 
estimated  at  10.2  minutes  (10).  This 
test  continued  to  be  conducted  under 
reduced  load  until  14.6  minutes,  when 
fire  exposure  was  terminated  without 
collapse  occurring.  The  predicted 
time-to-failure  falls  within  this  10-  to 
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15-  minute  range.  This  result  is  most 
promising  for  future  use  of  the  model. 

An  interesting  examination  is  how 
time-to-failure  is  altered  for  the  same 
truss  by  reducing  the  applied  load. 
The  failure-model  equation  predicts 
times-to-failure  as  a  function  of  ap¬ 
plied  load  as  shown  in  figure  5.  For 
reduction  in  load  to  50  percent  of  full 
design,  it  is  seen  5  minutes  is  added 
to  the  predicted  time  under  full 
design  load.  If  there  were  no  load  on 
the  floor  assembly  except  the  dead 
weight  (4.9  lb/ft2)  of  the  assembly 
itse'f,  a  failure  time  of  21.1  minutes 
could  result.  Hence,  failure  times 
greater  than  this  are  theoretically  im¬ 
possible  for  this  truss  design. 

Estimating  Floor- 
Assembly  Safety 

In  the  preceding  paragraphs  the 
authors  developed  models  for  predict¬ 
ing  the  time-to-failure  of  two  floor- 
assembly  types,  and  have  given  an 
accepted  model  to  predict  severity  of 
fire  exposure.  They  are  given  in  units 
of  time. 

Using  the  standard  reliability  nota¬ 
tion,  the  assembly  time-to-failure  and 
the  time  duration  of  fire  exposure,  tf_ 
and  t^  are  denoted  as  follows: 


R  =  tf 

(14) 

S  =  fj 

(15) 

variables  W,  Fj  A,  and  H.  The  real 
problem  lies  in  the  fire  severity  of 
E-119  not  being  representative  of  the 
fire  severities  associated  with  actual 
fire  situations. 

It  is  thus  appropriate  in  this  analysis 
to  treat  R^and  S^as  independent  ran¬ 
dom  variables. 

Using  the  approach  of  Zahn  (42), 

Re  I  =  Pr(R/S  >  1)  (17) 

and  taking  the  logarithm  of  the 
arguments 

Rel  =  Pr(ln(R/S)  >  0)  (18) 

By  making  the  following  definition 

J  =  ln(R/S)  (19) 

there  results 

Rel  =  Pr(J  >  0)  =  1  -  Fj  (0)  (20) 

where  FjJs  the  cumulative  density 
function  of  the  variable 
Using  first-order,  second-moment 
approximations,  the  mean  and  the 
variance  are  given  by 

MjSIn  -^3-  (21) 

°J3QR+QS  (22) 


Z  =  h(X„  X„. .  ,Xn)  (26) 
is  given  by 

E(Z)  =  h[E(X,),E(X)2), .  .  ,,E(Zn))  (27) 

Performing  this  operation  on  equation 
(6)  for  the  conventional  joist-floor 
assembly  and  replacing  the  expected 
values,  E,  of  the  component  variables 
by  their  statistical  estimates  denoted 
by  a  superscript  bar,  the  result  is 

^  2Cd(d  +  b)  +  6WKy/B 

MR  =  - - - — - 

V{2  Cd(d  +  b)  +  6MKr/B>J 

-  4C2  (b  +  4d)(bd;  -  6M/B) 
2C2(b  +  4d) 

(28) 


a  ~  V'b?  -  4  a,cT  (29) 

2a, 

where 

a^=  4CJ  (b  +  2d)  +  2PC  y  KfT  (29a) 

b,  ~  2PC/T-2dC(d  +  2b)  (29b) 

-  Pd  y  K/T  -  6M  y  K/B 

c,  =  bd2  -  Pd/T  -  6M/B  (29c) 


With  R^and  ^defined,  the  fire 
reliability  (Rel)  of  an  exposed  joist  or 
floor  truss  (as  a  part  of  a  floor),  given 
a  fully  developed  fire,  is  given  by 

Rel  =  Pr(R  >  S)  (16) 

which  reads  the  probability  (Pr)  that 
the  fire  resistance  is  greater  than  that 
of  the  fire  load.  It  is  convenient  for 
calculations  R^and  S^to  be  indepen¬ 
dent.  In  other  words,  the  fire  load  and 
associated  parameters  cannot  be  cor¬ 
related  to  the  members'  resistance, 
char  rate,  and  other  factors.  Con¬ 
versely,  structural  load,  member 
resistance,  and  char  rate  cannot  be 
influenced  by  fire  load.  Unfortunately, 
it  is  known  that  char  rate  is  cor¬ 
related  with  the  variables  of  equation 
(1)  that  have  been  discussed  by 
Schaffer  (34).  However,  in  practice, 
the  joist  or  component  will  be  ex¬ 
posed  to  a  standard  fire  condition 
such  as  ASTM  E-119  in  which  char 
rate  will  not  be  influenced  by 


Standardizing  J, 


°J 

and 

Rel  =  Pr  |  A>  -  —  j 

If  ttio  distribution  of  X  is  similar  in  all 
applications,  then  the  variable 

In 

ft  =  (25) 

°j  VQR+  SS 

is  a  consistent  measure  of  fire  safety. 
pjs  normally  called  the  safety  index. 

The  next  step  is  to  estimate  the 
means  and  the  variances  of  the 
resistance  and  the  load,  pg  and  ng,  o£) 
and  oi$  A  first-order  approximation  of 
the  mean,  E,  of  a  function,  Z,  where 


(23) 


(24) 


Again,  ^Js  treated  as  a  random 
variable  and  B^and  T^are  treated  as 
constants  equal  to  the  average  bend¬ 
ing  and  tensile  strength  of  the  lumber 
grade. 

The  same  operation  is  performed 
on  the  fire-severity  equation,  which 
results  in  _ 

W  Ap 

n^?nr  (30» 

The  first-order,  second-moment  ap¬ 
proximation  of  the  variance  of  Z 
defined  by  equation  (26)  is  given  by 


(31 1 


provided  the  Xj's  are  uncorrelated. 
The  partials  are  evaluated  at  their 
respective  mean  values. 
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In  the  expression  for  the  load 
variable  ^  the  floor  area  Ap  and  win¬ 
dow  area  A yv  could  be  correlated 
because  it  would  be  expected  Ayv  be 
some  way  related  to  the  perimeter 
that  involves  the  same  variables  as 
the  floor  area.  However,  because  data 
are  not  available  to  substantiate  this 
correlation,  a  zero  correlation  will  be 
assumed.  In  light-frame  construction, 
the  other  pairs  of  variables  lack  an 
obvious  cause  for  correlation. 

Without  showing  the  computations 

°l> =  SMCW+ 

+  °AW  +  Oj*/4>  (32, 


Qli  =  Qfa  +  Q\f 

+  °AW  +  2^4  (33) 

where  the  Q's  are  the  respective  COV. 

In  the  equation  for  resistance,  or 
time-to-failure,  of  the  floor  joist,  no 
obvious  confounding  correlations  ap¬ 
pear.  Using  equation  (31)  to  estimate 
the  variance  of  equation  (6),  three  par¬ 
tial  derivatives  must  be  calculated. 
These  derivatives  are  lengthy  so  they 
will  only  be  substituted  into  equation 
(31)  symbolically,  which  results  in 


Again,  the  partial  derivatives  are 
evaluated  at  the  mean  values  of  the 
component  variables. 

For  the  floor  truss,  component 
variables  M  and  F^are  perfectly  cor¬ 
related  with  a  correlation  of  +  1. 
Therefore,  equation  (31)  must  be 
altered  to  include  correlations  as 


i<  j 


where  ejj  is  the  correlation  coefficient 
between  variables  Xj  and  Xj.  By 
assuming  zero  correlation  tor  the 
other  variables,  as  for  the  floor  joist, 
application  of  equation  (35)  to  equa¬ 
tion  (29)  yields 


The  variance  of  the  char  rate,  oc2, 
has  been  estimated  by  reanalyzing 
data  previously  developed  (32).  The 
mean  charring  rate,  C,  and  estimated 
variance,  oc!,  under  ASTM  E-119  fire 
exposure  for  coast  Douglas-fir  and 
southern  pine  are 

coast  Douglas-fir 

C  =  0.0245  in. /min,  oc2  =  5.56  x 
10'6 

southern  pine 

C  =  0.0299  in./min,  oc2  =  3.80  x 
10'6 

The  variance  in  the  strength- 
reduction  factor,  y^is  unavailable. 
Substantial  information  is  available 
on  the  variation  of  the  applied  load 
that  defines  the  variation  of  N4  and  P. 
Eventually,  all  the  variances  denoted 
by  o2  will  be  replaced  by  statistical 
estimates  from  the  data  available  or 
from  the  results  of  present  research. 
Finally,  the  COV  of  Rcan  be  ob¬ 
tained,  and  is  given  by 


Qr  =  —  (37) 

MR 

As  statistical  data  become  available 
for  the  various  parameters,  the  com¬ 
ponents  of  the  safety  index  equation 
(25)  may  be  defined,  and  comparisons 
of  assembly  safety  accomplished. 


Discussion 

Probability  of  Failure 

Knowing  the  distribution  of  j^of 
equation  (23),  the  probability  of  struc¬ 
tural  failure  of  an  exposed  floor 
assembly,  given  the  occurrence  of  a 
fully  developed  fire,  can  be 
calculated.  In  general,  the  distribution 
of  z_is  not  known,  and  some  assump¬ 
tion  about  it  must  be  made.  Common 
practice  is  to  assume  ^follows  a  nor¬ 
mal  distribution.  In  this  analysis,  the 
assumption  of  normality  will  be  used 
realizing  the  probability  estimated  will 
be  in  the  neighborhood  of  0.1;  thus 
deviations  from  normality  from  one 


application  to  the  next  will  not  be 
amplified.  Under  these  assumptions, 
the  probability  of  failure,  Pf,  is  given 
by 

P,  =  <*>(-/?)  =  1  -  <t>(/J)  (38) 

where  <t>  is  the  cumulative  area  under 
the  standard  normal  curve. 

Code  Calibration 

In  recent  years  the  use  of  engineer¬ 
ing  components  has  increased 
dramatically.  Often  these  components 
are  fabricated  with  manmade 
materials,  and  the  variability  of  the 
mechanical  properties  of  these 
materials  is  substantially  less  than 
the  variability  of  a  natural  material, 
such  as  wood.  The  shortcoming  of 
the  present  "fire-rating”  system  is  it 
does  not  account  for  the  variability  of 
the  component  response,  although  it 
may  account  for  the  average 
response  of  a  component  to  fire. 

Quite  simply,  the  present  system  of 
fire  rating  allows  using  two  different 
components  with  an  equal  "fire 
rating"  of,  say,  1  hour,  but  at  the 
same  time  the  components  have  une¬ 
qual  safety  levels. 

The  situation  of  two  components 
with  unequal  safety  levels  can  best 
be  illustrated  by  a  hypothetical  exam¬ 
ple  using  equation  (25).  The  example 
involves  calculating  the  safety  index 
for  two  different  components,  the  on¬ 
ly  difference  being  the  variability  of 
the  time-to-failure  or  resistance  of  the 
component. 

For  the  example,  assume  that  the 
following  data  in  table  5  apply  to  two 
component  types  A  and  B.  Each  col¬ 
umn  is  identical  except  for  the  col¬ 
umn  indicating  the  COV  of  fire 
resistance.  On  applying  equation  (25) 
to  each  type,  the  resulting  safety  in¬ 
dices  are  jl^  =  0.98  and  /J3  =  1.24, 
which  indicates  an  unequal  level  of 
safety.  Converting  these  |Ts  by  equa¬ 
tion  (38)  to  probabilities,  the  results 
are  0.163  for  component  A  and  0.107 
for  component  B.  It  may  be  hastily 
argued  that  no  difference  exists  be¬ 
tween  a  16.3  percent  and  a  10.7  per¬ 
cent  failure  rate,  but  on  closer  ex¬ 
amination,  calculations  show  compo¬ 
nent  B  could  have  a  fire-endurance 
time  of  51.9  minutes  and  have  the 
same  relative  safety  as  component  A 
at  60  minutes! 

The  hypothetical  situation  il¬ 
lustrates  just  one  possibility  of  how 
the  safety  index  could  be  used  to  ob¬ 
tain  equal  safety  and  give  a  "fair 


Table  5.— Level  of  safety  of  components,  A  and  B,  with  same  average  fire-endurance  time 
(60  min)  but  different  variabilities  (COV  of  0.25  versus  0.5)' 


Component 

CR 

2S 

Min 

Min 

A 

60 

0.5 

30 

0.5 

B 

60 

.25 

30 

.5 

'  Each  component  is  exposed  to  an  identically  distributed  tire  load,  as  illustrated  by  last  two  col¬ 
umns.  result  is  an  unequal  level  of  safety,  as  calculated  by  the  safety  index  equation  (25). 


shake”  to  new  materials  and  new 
components. 

In  the  future  it  may  be  possible  to 
identify  at  what  point  fire  load,  hence 
fire  duration,  have  different  expected 
values  end  different  levels  of  variabili¬ 
ty.  In  this  type  of  situation,  a  lower  or 
a  higher  fire-rated  component  may  be 
needed. 

Summary  and 
Conclusions 

1.  Time-to-structural-failure  predic¬ 
tion  models,  based  on  the  residual 
load-carrying  capacity  of  fire-exposed 
floor  elements,  are  given  for  two  un¬ 
protected  light-frame  wood  floor 
assemblies.  A  reduction  in  strength 
factor,  a^was  calculated  from  the 
limited  fire-exposure  test  data  accord¬ 
ing  to  two  equations: 


Joist  floor 

“Ml  =  oB 
Ktf.C) 

Floor  truss 

fb  x  ft 

—  +  —  =  a 

B  T 

The  thermal  reduction  factor  is  fur¬ 
ther  a  function  of  fire-exposure  time 
and  the  geometry  of  small  cross  sec¬ 
tions: 

Joist  floor 

a  =  1/(1  +  yKtf) 

where 

y  =  0.170  (in./min) 

Floor  truss 

a  =  i/(i  +  yKtf) 


where  for 

Douglas-fir 

y  =  0.113  in./min 

Southern  pine 

y  =  0.0839  in./min 

2.  A  comparison  of  predicted 
times-to-failure  versus  those  actually 
observed  for  four  unprotected  joist- 
floor  assemblies  results  in  predicted 
times  consistently  less  than  those 
observed.  The  difference  is  attributed 
to  three  factors: 

—The  model  has  parameters  quan¬ 
tified  on  the  basis  of  fire-endurance 
tests  of  paired  joists  with 
negligible  floor  sheathing. 

—The  actual  floors  consist  of  many 
joists  with  load  sharing  likely. 

—The  actual  floors  have  floor 
sheathing  that  contributes  to  in¬ 
creased  load-carrying  capacity. 

Predicted  times-to-failure  for  two  of 
the  floor  assemblies  were  similar  to 
those  of  observed  times-to-failure  of 
the  first  joist  (not  total  floor  assembly 
failure). 

3.  The  time-to-failure  model  with  in¬ 
put  parameters  for  a  southern  pine 
floor-truss  assembly  that  had  been 
fire-endurance  tested  resulted  in  a 
predicted  time-to-failure  of  11.2 
minutes.  The  fire  test  had  been  con¬ 
cluded  at  10.2  minutes  because  of  ex¬ 
cessive  deflection  of  the  floor 
assembly,  without  evidence  of  col¬ 
lapse. 

4.  The  time-to-failure  model  for  a 
loaded  and  a  fire-exposed  floor-truss 
assembly  was  employed  to  examine 
the  influence  of  various  floor  loads  on 
time-to-failure.  At  full-design  load, 
failure  was  predicted  to  be  11.2 
minutes;  but  with  only  dead  load,  the 
period  was  extended  to  21.1  minutes. 
This  indicates  the  sensitivity  of  this 
type  of  fire-exposed  assembly  to  the 


load  applied  during  test.  In  this 
analysis  it  was  critically  assumed 
that  failure  of  the  lower  chord,  rather 
than  connectors,  would  result  in  col¬ 
lapse  of  the  assembly. 

5.  The  procedure  to  use  for 
calculating  safety  of  unprotected 
light-frame  floor  assemblies  is  given. 
An  example  is  provided  to 
demonstrate  how  variability  in 
assemblies  can  affect  the  com¬ 
parative  safety  of  assemblies. 

6.  The  predictive  capabilities  for 
both  of  the  proposed  assembly 
models  require  further  fire-exposure 
experiments  for  independent  valida¬ 
tion  and  parameter  refinement. 


Nomenclature 


Ap  Floor  area 

A yy  Window  or  opening  area 

b  Beam  breadth 

B  Modulus  of  rupture  at  room 

temperature 
C  Char  rate 

d  Beam  depth 

F  Allowable  stress 

H  Window  or  opening  height 

I  Moment  of  inertia 

K  Ratio  of  perimeter  to  area  of 

cross  section. 

M  Applied  moment 

MOE  Modulus  of  elasticity 
P  Axial  load 

R  Member  or  structural  resistance 
S  Applied  “load” 

T  Ultimate  tensile  stress 

to  Fire  duration 

W  Fuel  load  density 

Y  Distance  from  beam  centroid  to 

outer  fiber 

a  Ratio  of  high-temperature  to 

normal-temperature  strength 
p  Safety  index 

y  Thermal  degrade  factor 

fj  Statistical  mean 

f  Correlation  coefficient 

o  Standard  deviation 

0  Accumulative  density  function 
for  standard  normal  distribution 
Q  Coefficient  of  variation 

Subscripts 

b  Bending 

d  Duration 

f  Failure 

R  Resistance 

S  Load 

t  Tension 


12 


Literature  Cited 


1.  Abrams,  M.  S. 

1978.  Behavior  of  inorganic  materials  in  fire.  ASTM  E-5  Symposium  on 
design  of  buildings  for  fire  safety.  Boston  Park  Plaza,  Boston,  Mass. 

2.  American  Iron  and  Steel  Institute. 

1978.  Designing  fire  protection  for  steel  columns.  American  Iron  and  Steel 
Inst.,  Washington,  D.C. 

3.  American  Society  for  Testing  and  Materials. 

1977.  Standard  methods  of  fire  tests  of  building  construction  and 
materials  (E-119).  Annual  Book  of  ASTM  Standards,  Part  18.  ASTM, 
Philadelphia,  Pa.  p.  739-757. 

4.  Ang,  A.  H.S.,  and  others. 

1972.  Structural  safety— A  literature  review.  J.  Structural  Division,  Proc. 
ASCE  98  (ST4):845-884. 

5.  Burros,  R.  H. 

1975.  Probability  of  failure  of  building  from  fire.  J.  Structural  Division. 
Proc.  ASCE  101  (ST9):1947-1960. 

6.  Coward,  S.K.D. 

1975.  A  simulation  method  for  estimating  the  distribution  of  fire  severities 
in  office  rooms.  Build.  Res.  Establ.  Current  Pap.  31/75.  Fire  Res.  Stn., 
Borehamwood  Hertfordshire,  England,  WD6  2BL. 

7.  Doyle,  D.  V.,  and  L.  J.  Markwardt. 

1966.  Properties  of  southern  pine  in  relation  to  strength  grading  of  dimen¬ 
sion  lumber.  U.S.  Dep.  Agric.,  For.  Serv.  Res.  Pap.  FPL  64.  Forest  Products 
Laboratory,  Madison,  Wis. 

8.  Doyle,  D.  V.,  and  L.  J.  Markwardt. 

1967.  Tension  parallel-to-grain  properties  of  southern  pine  dimension 
lumber.  U.S.  Dep.  Agric.,  For.  Serv.  Res.  Pap.  FPL  84.  Forest  Products 
Laboratory,  Madison,  Wis. 

9.  Ellingwood,  B.,  and  J.  R.  Shaver. 

1977.  Reliability  of  RC  beams  subjected  to  fire.  J.  Structural  Division, 

Proc.  ASCE  103  (ST5):1047-1059. 

10.  Factory  Mutual  Research. 

1977.  Fire  endurance  test  of  floor-ceiling  assembly.  Wood  trusses  with 
plywood  floor  design.  ASTM  E-1 19-76,  FC-250.  Factory  Mutual  Res. 

11.  Gross,  D. 

1977.  Measurements  of  fire  loads  and  calculations  of  fire  severity.  Wood 
and  Fiber  9(1):72-85. 

12.  Gustaferro,  A.  H.,  and  D.  P.  Jenny. 

1978.  Alternative  to  fire  testing— Design  of  concrete  structures  for  fire 
endurance.  Southern  Building.  June-July:24-28. 

13.  Gypsum  Association. 

1978.  Fire  resistance  design  manual.  Gypsum  Assoc.,  Evanston,  III. 

14.  Hoyle,  R.  J. 

1976.  Digest  of  tension  parallel-to-grain  strength  of  Inland  North  Douglas- 
fir.  Res.  Rep.  76/57-22.  Coll,  of  Eng.  Res.  Div.,  Washington  State  Unlv., 
Pullman,  Wash. 

15.  Hoyle,  R.  J. 

1977.  Review  of  world  literature  on  characteristic  distribution  of  proper¬ 
ties  of  2-inch  softwood  dimension  lumber.  Res.  Rep.  77/57-11.  Coll,  of  Eng 
Res.  Div.,  Washington  State  Univ.,  Pullman,  Wash. 

16.  Hoyle,  R.  J.,  and  T.  M.  Maloney. 

1976.  Bending  strength  tests  of  visually  graded  2-inch  dimension  lumber 
from  western  Canada.  Res.  Rep.  No.  76/57-36.  Coll,  of  Eng.,  Res.  Div., 
Washington  State  Univ.,  Pullman,  Wash. 


I 


17.  Hoyle,  R.  J.,  and  T.  M.  Maloney. 

1976.  Tension  strength  and  bending  elastic  modulus  of  truss  joist 
machine  stress  rated  2  by  4  lumber.  Res.  Rep.  76/57-44.  Coll,  of  Eng.  Res. 
Div.,  Washington  State  Univ.,  Pullman,  Wash. 

18.  Kameda,  H.,  and  T.  Koike. 

1975.  Reliability  theory  of  deteriorating  structures.  J.  Structural  Division, 
Proc.  ASCE  101  (ST1):295-310. 

19.  Knudsen,  R.  M.,  and  A.  P.  Schniewind. 

1975.  Performance  of  structural  wood  members  exposed  to  fire.  Forest 
Products  J.  25(2):23-32. 

20.  Lawson,  D.  I. 

1952.  The  fire  endurance  of  timber  beams  and  floors.  Structural  Engineer, 
30(3):27-33. 

21.  Lie,  T.  T. 

1972.  Fire  and  buildings.  Applied  Science  Publishers,  Ltd.  Ripple  Road, 
Barking,  Essex,  England.  276  p. 

22.  Lie,  T.  T. 

1972.  Optimum  fire  resistance  of  structures.  J.  Structural  Division,  Proc. 
ASCE  98(ST1):2 15-232. 

23.  Lie,  T.  T. 

1974.  Characteristic  temperature  curves  for  various  fire  severities.  Fire 
Technology  10(4):3 15-326. 

24.  National  Forest  Products  Association. 

1974.  ASTM  E-119  fire  endurance  test— 2-  by  10-inch  wood  joist  floor 
assembly.  Design  FC  209.  Nat.  For.  Prod.  Assoc.,  Washington,  D.C. 

25.  National  Forest  Products  Association. 

1974.  ASTM  E-119  fire  endurance  test— 2-  by  10-inch  wood  joist  floor 
assembly.  Design  FC  212.  Nat.  For.  Prod.  Assoc.,  Washington,  D.C. 

26.  National  Forest  Products  Association. 

1974.  ASTM  E-119  fire  endurance  test— 2-  by  8-inch  wood  joist  floor 
assembly.  Design  FC  213.  Nat.  For.  Prod.  Assoc.,  Washington,  D.C. 

27.  National  Forest  Products  Association. 

1974.  ASTM  E-119  fire  endurance  test— 2-  by  8-inch  wood  joist  floor 
assembly.  Design  FC  216.  Nat.  For.  Prod.  Assoc.,  Washington,  D.C. 

28.  National  Forest  Products  Association. 

1977.  Design  values  for  wood  construction— a  supplement  to  the  1977 
edition  of  national  design  specification  for  wood  construction. 
Washington,  D.C.  20036.  20  p. 

29.  Raes,  H. 

1977.  The  influence  of  a  building's  construction  and  fire  load  on  the  inten¬ 
sity  and  duration  of  a  fire.  Fire  Prevention  Sci.  and  Technol.  16:4-16. 

30.  Robertson,  A.  F„  and  D.  Gross. 

1970.  Fire  load,  fire  severity,  and  fire  endurance.  Special  Technical 
Publication  464.  American  Society  for  Testing  and  Materials,  Philadelphia, 
Pa.  p.  329. 

31.  Schaffer,  E.  L. 

1961.  The  effects  of  fire  on  selected  structural  timber  joints.  M.S.  thesis, 
Univ.  of  Wis.,  Madison,  Wis. 

32.  Schaffer,  E.  L. 

1966.  Charring  rate  of  selected  woods— transverse  to  grain.  U.S.  Dep. 
Agric.,  For.  Serv.  Res.  Pap.  FPL  69.  Forest  Products  Laboratory,  Madison, 
Wis. 

33.  Schaffer,  E.  L. 

1973.  Effect  of  pyrolytic  temperatures  on  the  longitudinal  strength  of  dry 
Douglas-fir.  ASTM  J.  Test  Evaluation  1(4):319-329. 


34.  Schaffer,  E.  L. 

1977.  State  of  structural  timber  fire  endurance.  Wood  and  Fiber 
9(2):  145-1 70. 

35.  Son,  B.  C. 

1973.  Fire  endurance  tests  on  unprotected  wood-floor  constructions  for 
single-family  residences.  Center  for  Build.  Tech.  Inst,  for  Applied  Technol., 
National  Bureau  of  Standards,  NBSIR  734263.  Washington,  D.C. 

36.  Suddarth,  S.  K. 

1972.  A  computerized  wood  engineering  system:  Purdue  plane  structures 
analyzer.  U.S.  Dep.  Agric.,  For.  Serv.  Res.  Pap.  FPL  168.  Forest  Products 
Laboratory,  Madison,  Wis. 

37.  Suddarth,  S.  K.,  F.  Woeste,  and  W.  Galligan. 

1978.  Differential  reliability:  Probabilistic  engineering  applied  to  wood 
members  in  bending/tension.  U.S.  Dep.  Agric.,  For.  Serv.  Res.  Pap.  FPL 
302.  Forest  Products  Laboratory,  Madison,  Wis. 

38.  Sunley,  J.  G. 

Design  concepts  for  fire-resisting  constructions.  Behavior  of  wood  pro¬ 
ducts  in  fire.  Pergamon  Press,  Oxford,  p.  95-101. 

39.  Thomas,  P.  H.,  and  C.  R.  Theobald. 

1977.  Part  2:  The  burning  rates  and  durations  of  fires.  Fire  Prevention  Sci. 
and  Technol.  17:15-16. 

40.  U.S.  Department  of  Commerce/National  Bureau  of  Standards. 

1978.  Building  research  translation— the  behavior  of  concrete  structures 
in  fire— a  method  for  prediction  by  calculation.  NBS  Tech.  Note 
710-7110.  U.S.  Government  Printing  Office,  Washington,  D.C. 

41.  U.S.  Department  of  Housing  and  Urban  Development. 

1973.  HUD  minimum  property  standards— one-  and  two-family  dwellings. 
U.S.  HUD,  Washington,  D.C. 

42.  Zahn,  J.  J. 

1977.  Reliability-based  design  procedures  for  wood  structures.  Forest  Pro¬ 
ducts  J.  27(3):21-28. 


to 

•H 

3 

to 

•» 

O  C  T3 

c  • 

4T  G)  rH 

O  '""v 

4J  E  P  «-l 

W  CN 

•H  C  O  CO 

-<r 

u-i  o  to  >  <u 

Jn'O  M 

o  o  -H 

tO  «  o 

a  to  e  u  * 

R  l 

c  K  y  c  w  e 

•  »J 

c  «H  C  o  c 

•>  *  fc- 

-H  C  p  V  CX-H 

XI  6-1  lx 

y  c3  «  o  v 

O  V 

B3  4)  t»  *0  y  r3 

oca; 

3  f-4  u  v  a  N 

3  C  -M 

*H  O  c  -H 

>1 

■a  o 

CO  G>  «3  -O  rH 

X. 

C  32 

S>  *  *5  P  O  -H 

o 

CJ  H  4J  rH  C  4-J 

-H  O  * 

H  £  SO  ^  3 

?3 

43  *3 

C  «G  H 

X, 

e  c  c 

C  E  C  C  <0 

o 

O  -h  us 

to  o  p  > 

I-H  yj 

«*-i  •  *H  *C  «H 

cz 

O  Ci  • 

o  v  io  y  y  y 

J 

u  «  > 

o  *o  to  o  u 

X  — 

to  c  £  cj 

w 

C  P 

XI  to  O  k  60  <M 

*-* 

c  *a  to 

rH  l>  ?  O  C  v 

f>  P 

P  VH  *H  *ri  t) 

i)  «J  • 

CO  h  C  C  C 

p-! 

*3  C 

GJ  C  <0  **-l  fc-i  O 

c  o 

U  <U  «H  tl  y 

Vl  QJ  ft, 

atfj  ^  cj 

t— * 

O  <0 

«u  o  B  *n  c  t-* 

■u  < 

*c  u  o  o  c 

10  X  O 

y  D.H’O  u  d 

HI  C  w 

o  c  yj 

■n  #3 

to  h  u  y  o  v-i 

o 

W  «  W 

w  (0  C  C  O 

l-l 

C  O  C  P  3  3. 

•-H 

G>  C  * 

0)  i-i  <u  e  4-i  E 

• 

a  c  a 

(n  c  >  E  is  n 

/. 

O  rH 

C  CTJ  <u  ^  ^ 

• 

t-4  <  CM 

l-i  _£L  to  (X  P  O 

P-i  *H 

p-  o  ex  *j  to 

• 

(U  >y  D  •H  H 

w 

E  O  to  n  « 

03 

*H 

3 

•3 

•> 

a;  c  T3 

c  • 

£  <U  rH 

o  ^ 

4->  g  oi  x 

•5  CN 

■H  >-i  O  tO 

+4  <r 

14-1  o  is  S  p 

CN 

o  a;  vt 

•o  «  o 

p.  k  E  y  • 

£  i 

c  <0  -u  O  Cl  c 

-  ►j 

O  H  y  4)  o 

03  -  {X 

n  >4  3  si  an 

*T3  Vi  U* 

y  to  k  o  -u 

O  v 

IS  p  P  T3  6i  IS 

o  o  a> 

d  h  p  d  a  n 

3  C  -U 

rH  O  C  *H 

•o  0 

«  3  #S  H 

C  3  2 

>  *  to  P  O  vt  i 

is  .e 

4i  H  il  H  o  y 

*ri  O  • 

H  £1  U1  ?  3 

ro 

»Q  tn 

C  «  H 

M 

ECO 

CO  B  c  Cl  o 

o 

O  *H  &. 

03  O  C  > 

fH  4-1 

Vi  •  *H  V4 

o  u  • 

O  <M  to  4J  4J  4J 

tJ  <0  > 

o  -a  «o  o  u 

£  t- 

to  o  fc  a» 

f/l 

c  <o 

U  to  O  M  60  6-< 

1)  T3  LO 

^4  (U  >  O  C  v 

-J 

>  c 

3  «H  V  *H  C 

— ' 

o  «  • 

to  4-4  C  C  p 

- 

t/3  c 

<y  C  «J  «H  C  O 

c  o 

y  i  <h  c  h 

IM  Si  3- 

aj?  >.  u 

_ _ , 

o  03 

P  O  f  rfl  C  4J 

4-*  < 

x:  c  c  o  c 

n  -o  c 

y  an  o  u  <4 

O  C  00 

c  C  H 

v  a> 

91  H  U  y  4)  y 

U  C2 

h  <0  c  c»  o 

c 

p  y  c  p  3  a 

— 

a;  c  • 

P  1-4  P  g  4-1  E 

o.«  a 

to  p  >  p  n  -H 

/. 

P  IS  P  -4  X 

6i  <  in 

y  £  fi  cuj  o 

a  _. 

a-  *-» 

pc  u  C  y  to 

^  Vt  3  "H  H 

00 

E  O  tfl  H  fl 

(0 

*H 

& 

(0 

» 

<13 

e 

T3 

a 

• 

*e 

0) 

rH 

o 

e 

GJ 

c< 

to 

CN 

Ci 

o 

03 

Ti 

<r 

V 

O 

<0 

4) 

>>  T3 

CN 

o 

0) 

*H 

to 

o 

a 

CO 

£ 

4-1  • 

x: 

l 

c 

0) 

4-» 

o 

c*  c 

hJ 

o 

rH 

Cl 

4>  O 

to  - 

p-i 

•H 

Cl 

2 

V 

(X  «H 

»o  v 

fy 

<U 

10 

03 

O  4-» 

O  V 

<0 

o 

V 

•o 

Cl  CO 

o  o 

a) 

3  rH 

Cl 

0) 

CX  N 

>> 

>  c 

4-> 

rH 

u 

c 

Cl 

*o 

o 

cO 

0) 

CO 

TJ  «H 

o 

x: 

> 

•• 

CO 

(U 

O  t4 

4J 

«S6 

a> 

rH 

0) 

*H 

2  ** 

CO 

• 

rH 

& 

to 

?  P 

Cl 

43 

to 

p 

«0 

H 

O 

e  (3 

€ 

cO 

P 

J2 

c  tu 

A 

O  -H 

p£ 

CO 

o 

<u  > 

CO 

•H  4J 

V 

• 

-H 

*C  *H 

o  c 

• 

o 

V 

03 

4J 

4J  4J 

o  «o 

> 

o 

*t3 

tO 

O  U 

(0 

X 

u 

CO 

O 

£ 

o 

4-1 

o 

c 

<u 

4J 

03 

O 

Ci 

oy  V 

a>  *o 

to 

»H 

0) 

> 

o 

C  v 

1 

>  c 

0 

*H 

V 

*h  tu 

<U  <0 

• 

(0 

4J 

C 

c 

C 

tn 

u 

tu 

C 

cO 

*H 

c  o 

Ci 

a. 

a 

o 

c 

<u 

-H 

to  4-1 

V  <u 

p-i 

P-  40 

o 

tJ 

O  03 

QJ 

o 

E 

43 

C  4J 

C/3 

4J 

< 

42 

Cl 

o 

o  c 

<D 

CO  'O 

p 

cy  *h 

•o 

o  «0 

C< 

<u  c 

to 

O 

a; 

4-1 

o 

«H  0) 

p 

to 

•H 

CJ 

u 

a>  c 

Cc 

4-1  03 

s-/ 

O 

CO 

C 

c»  o 

C* 

C 

a 

c 

a> 

3  CX 

• 

0)  c 

w 

0) 

a; 

e 

c»  £ 

<Tj 

cc  cO 

p 

(0 

c 

3> 

S 

C0  *r4 

O  f-i 

a) 

to 

CL> 

rH 

U 

C»  <  CN 

c 

42 

03 

P 

to  o 

Pw 

1 — 1 

PC 

O 

P 

4J  03 

0) 

V 

D 

•H  t— 1 

PO 

£ 

o 

CO 

rH  CO 

CO 

•H 

03 

•« 

(0 

c 

*T3 

C  • 

42 

(0 

rH 

O  /-N 

4-> 

£ 

GJ 

Ci 

03  CN 

-H 

C 

O 

CO 

•h  <r 

V 

O 

CO 

£ 

G) 

^  n3  CN 

0 

to 

•H 

43  50 

cx 

CO 

6 

4-1  • 

X  1 

c 

03 

4J 

O 

Cl  <2 

-  rJ 

0 

rH 

Ci 

<0  O 

(0  •>  p 

•H 

Ci 

0 

V 

CX  *H 

•0  V  fcy 

4J 

«o 

CO 

O  4H 

O  VI 

<0 

to 

<0 

T3 

c»  <0 

0  0  4) 

0 

rH 

Cl 

<0 

CX  N 

^  c  y 

rH 

O 

c 

•H 

-0  0 

*0 

G) 

cO 

-O  rH 

Cl 

C  32 

> 

» 

03 

CJ 

O  *H 

o 

CO  42 

*w  0  • 

a> 

rH 

GJ 

iH 

6  4J 

■*-» 

iH 

S 

W3 

?  3 

cO 

43  03 

e 

CO 

Cl 

B  G  to 

to 

£ 

<2 

Cl  0) 

o 

O  *H  P2 

03 

O 

<0  > 

43 

rH  4H 

V 

• 

42  *H 

cd 

O  Ci  • 

0 

V 

CO 

4J 

4-1  4-> 

,J 

CJ  to  > 

O 

T> 

«0 

0  u 

X  c 

to 

n 

f 

<0 

CO 

{2  <0 

4-1 

to 

P 

Ci 

0C  V4 

<0  T3  CO 

rH 

0) 

5 

O 

C  V 

o 

>  c 

3 

-H 

V 

•H  GJ 

0)  «3  • 

03 

4-1 

c 

(2 

C 

E 

(0  Ci 

<0 

Cl 

CO 

■H 

C  0 

c  c 

Cl 

0 

•H 

CL)  4-i 

Ci 

v  <0  u- 

a  43 

>> 

U 

a. 

0  to 

<0 

0 

£ 

43 

C  4-1 

H  < 

42 

Ci 

O 

0  c 

to 

Cl 

to  TJ  P 

4-1 

CX  rH 

T3 

u  to 

<0  C  to 

O 

Qi 

4-1 

H  4)  D 

CO 

rH 

CJ 

4-1 

<0  Cl 

o 

y  co  ^ 

4J 

«0 

C 

Cl  0 

lu 

Ci 

c 

U 

c 

GJ 

3  9*  1 

G>  C  • 

a> 

•H 

0/ 

E 

H  6  ! 

• 

cx  to  CX 

03 

C 

> 

(0 

CO  -H  1 

C/3 

O  *H 

G> 

fl 

<0 

rH 

C 

Cl  <  CN 

U 

s: 

03 

CX 

GJ  O 

a-  rH 

CX 

u 

cx 

4J  0) 

■ 

0) 

Vi 

3 

■H  rH 

00 

6 

O 

0) 

rH  id 

